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Abstract: The radio technique for the detection of cosmic particles has seen a major revival in recent 
years. New and planned experiments in the lab and the field, such as GLUE, Anita, LUNASKA, Co- 
dalema, LOPES as well as sophisticated Monte Carlo experiments have produced a wealth of new in- 
formation and I review here briefly some of the main results with the main focus on air showers. Radio 
emission of ultra-high energy cosmic particles offers a number of interesting advantages. Since radio 
waves suffer no attenuation, radio measurements allow the detection of very distant or highly inclined 
showers, can be used day and night, and provide a bolometric measure of the leptonic shower component. 
The LOPES experiment has detected the radio emission from cosmic rays, confirmed the geosynchrotron 
effect for extensive air showers, and provided a good calibration fomula to convert the radio signal into 
primary particle energy. Moreover, Monte Carlo simulations suggest that also the shower maximum and 
the particle composition can be measured. Future steps will be the installation of radio antennas at the 
Auger experiment to measure the composition of ultra-high energy cosmic rays and the usage of the 
LOFAR radio telescope (and later the SKA) as a cosmic ray detector. Here an intriguing additional ap- 
plication is the search for low-frequency radio emission from neutrinos and cosmic rays interacting with 



the lunar regolith. This promises the best detection limits for particles above 10 
go significantly beyond current ground-based detectors. 



eV and allows one to 



Introduction 

Radio astronomy has always been closely con- 
nected to cosmic ray physics. Already the very 
first cosmic radio emission detected by Carl Jansky 
in 1932 originated from cosmic ray constituents 
in the MiUcy Way. We now know that at low ra- 
dio frequencies the diffuse Galactic radio emis- 
sion is mainly produced through synchrotron radia- 
tion of relativistic electrons. They are propagating 
through the interstellar medium and the Galactic 
magnetic field and were most likely accelerated in 
supernova explosions. 

Also, the brightest radio sources discovered there- 
after, like quasars and radio galaxies (active galac- 
tic nuclei) and supernova remnants, are today the 
main suspects for the origin of cosmic rays. So, 
without the advent of radio telescopes, we proba- 
bly would not know much about the non-thermal 
universe today. 



It is therefore no surprise that radio antennas were 
early on also considered for directly detecting cos- 
mic ray air showers. In fact the huge Lovell radio 
telescope in Jodrell Bank was initially built in or- 
der to detect cosmic ray radar reflection [14, 27] — 
it did not succeed but detected the radar reflection 
of Sputnik instead and made history. 

Radio detection of air showers has a number of 
advantages: the detector material itself, a simple 
wire, is cheap, radio emission is not absorbed in the 
atmosphere and can thus see the entire shower, and 
interferometric techniques should allow relatively 
precise localization. But does this work in practice 
and how does the radio signal actually look like? 

I will here mainly summarize some of the main re- 
cent results and not recall the entire history of this 
field. Here one can point to the well-known review 
by Allan from 1971 [1] and a brief summary of the 
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early results given by Falcke & Gorham [20]. Main 
points were the prediction of radio Cherenkov ra- 
diation by Askaryan [6, 7] and the discovery of CR 
related radio pulses through Jelley et al. in 1965 
[39] (a nice historical recount of the discovery was 
given by Trevor Weekes [67]). 

Despite many experimental problems at the time, 
quite a number of basic properties of air shower 
were established within a decade culminating in 
the empirical "Allan formula" [1]. 

A long hiatus of this field began in the 1970's, 
as witnessed by a quote from Alan Watson in his 
1975 ICRC rapporteur talk in Munich, where he 
observed that "Apart from work at 2 MHz which 
is planned for Yakutsk, it is clear that experimen- 
tal work on radio signals has been terminated else- 
where." 

Occasional attempts with single or few radio anten- 
nas at EAS-TOP and CASA/MIA [26] did not lead 
to further radio detections in the 1990's, making 
some colleagues (unjustifiably) even doubt the re- 
ality of the earlier results in private conversations. 

Only in recent years, the technique has seen an 
astounding revival. A good overview is probably 
found in references [59, 20, 21, 52]. 

Scientifically, this started with attempts to detect 
radio emission from neutrinos hitting the moon by 
Hankins, Ekers, & O' Sullivan [29] and Gorham 
et al. [25]. For air showers, the realization that 
the emission can be understood as geosynchrotron 
emission by Falcke & Gorham [20] and Huege & 
Falcke [32] also inspired new efforts. 

Technologically, the revival is certainly due to 
high-dynamic range digital radio receivers and 
post-processing capabilities that are now available. 
There is also a general revival in low-frequency 
radio astronomy as seen in a number of projects 
such as LOFAR[22], MWA[15], LWA[44] and 
GMRT[64]. 

In the following we will give a summary of some 
of the results in this field, with particular empha- 
sis on radio emission from air showers and results 
obtained with LOPES. 



Theory 

After the experimental realization in the late 1 960's 
that the Earth magnetic field is a factor in radio 
air shower emission, early theoretical modeling by 
Kahn & Lerche considered the Lorentz boosted lat- 
eral current induced by the geomagnetic field [41], 
an approach that has been revisited very recently 
by Werner & Scholten [68, 63]. 

A different approach was presented by Falcke, 
Gorham, and Huege [20, 32] where the radio emis- 
sion was explained in terms of "geosynchrotron 
emission". This approach takes an important ex- 
tra factor into account, namely the curvature of 
the trajectories of the individual electron/positron 
pairs in the geomagnetic field. The fact that the 
emission region is smaller than a wavelength (and 
optically thin) allows a relatively simple coherent 
addition of the radio waves. This "single-particle 
approach" makes it straightforward to combine 
the radio emission with Monte Carlo calculations 
[33, 36]. The overall level of the geosynchrotron 
component seems to be sufficient to explain the 
bulk of the observed radio emission [34]. After 
all, geosynchrotron subsumes most of the "lateral 
current interpretation". 

However, also in the geosynchrotron picture a 
couple of extra effects still need to be taken 
into account, such as the current induced by 
the change in charges through creation, annihi- 
lation and recombination[63, 49]. The recom- 
bination of electrons may also lead to an addi- 
tional Bremsstrahlung component [49]. Also, the 
non-zero refractive index of air and the original 
Askaryan effect through Cherenkov emission from 
the charge excess [18, 50] will play a role at some 
level. The static Coulomb contribution should also 
be looked at [50] as well as optical depth effects 
that could become relevant for air showers around 
10^" eV (proposed in the context of radio radar ex- 
periments [23]). 

Hence, while major progress has been made, there 
is still room for improvement. Nonetheless, the 
predictive power of current Monte Carlo codes, if 
coupled with air shower simulations, is probably 
already quite significant. This requires knowledge 
of the lepton evolution in the showers and adequate 
shower libraries [47]. 



30th International Cosmic Ray Conference 



1e+11 



1e+10 



1e+09 



1e+08 




1000 10000 100000 

32 to 64 MHz amplitude at 275 m [nV m"^ 



1e+06 



1100 
1050 
1000 
950 
900 
850 
800 
750 
700 
650 
600 

















' *= 


)< 


_ 








> 


*\ 






_ 


- 






-H- 
+ 




X 






- 






^^ 


Jk^+ 










" 


+^ 


S»^ 


i * 










" 


_ 


^ 


^t^ 












_ 




ti 


^ 
















^^ 






iron 




m^" 


pV 






^r 




+ 


proton 


X 


< 


eV 




i 


1 




+ 
1 


gamma 

1 1 


)K 


10^° 

1 


eV 

1 


- 



10 



,/E7 



, at 32 to 64 IVIHz 



Figure 1: Integrated 32-64 MHz radio flux mea- 
sured at 275 m from REAS2 simulations as func- 
tion of energy for different primary particles and 
showers with 60° zentith angle. Shower-to-shower 
fluctuations are only 5%. Figure taken from refer- 
ence [35]. 



Figure 2: Location of shower maximum, X^ax, vs. 
lateral slope of radio emission from REAS2 simu- 
lations. The tight correlation suggests that perhaps 
Xmax should be measurable with radio antennas 
and showers with 60° zentith angle. Figure taken 
from reference [35]. 



Calculations with the REAS2 code and CORSIKA 
code by Huege et al. have recently shown some 
interesting results [35]: If measured at a charac- 
teristic radius of ~ 300 m from the shower core 
the radio signal is tightly correlated with the pri- 
mary particle energy (Fig. 1). Shower-to-shower 
fluctuations and different elemental composition of 
the primaries induce just 5% variations in the ra- 
dio flux. This is due to the fact that the radial ra- 
dio distribution on the ground pivots around a few 
hundred m for different X^ax, depending a bit on 
shower geometry. This also means that measur- 
ing the radial slope of the radio emission should 
give clues for the location of the shower maximum 
(Fig. 2) and together with the absolute radio flux al- 
low one to separate primaries of different elemen- 
tal composition. This tantalizing prediction nat- 
urally requires experimental confirmation but al- 
ready shows how important the theoretical work is. 

LOPES 

Quite a few experiments have been built or have 
been discussed in the last couple of years that em- 
ploy the radio detection method. We will try to 
briefly discuss them here in turn. 



A very productive experiment has been the LO- 
FAR PrototypE Station (LOPES), which made use 
of early prototype hardware developed for the LO- 
FAR radio telescope (see below) and helped to 
bring about the current renaissance in radio detec- 
tion techniques [19]. 

LOPES [30, 19] was a collaboration of radio as- 
tronomers involved in LOFAR and the groups 
involved in the KASCADE[45] and KASCADE 
Grande array [53]. The idea was to put a signifi- 
cant amount of radio antennas - allowing for inter- 
ferometric measurements - near a well-developed 
air shower array. This facilitates a cross corre- 
lation between conventional air shower measure- 
ments and radio observations. 

LOPES consisted initially of 10 single dipole an- 
tennas (LOPES 10) that were then expanded to 30 
antennas (LOPES30). In the last phase LOPES 
was again rearranged to have 20 antennas of which 
10 are in a dual polarization mode (Fig. 3). Polar- 
ization investigations are now underway (see Isar 
etal.[37]). 

The LOPES antennas digitize the incoming radio 
waves with a 12 bit A/D converter operating at 
80 MHz. An analog filter restricts the observ- 
able frequency range to 40-80 MHz, i.e. the second 
Nyquist zone. All antennas share a joint clock dis- 
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tribution and have a 6.7 second ring-buffer which 
is triggered and read-out roughly twice a minute by 
the KASCADE array (Fig. 3). 

LOPES itself is restricted to the dimensions of 
KASCADE (200 m), but with the help of KAS- 
CADE Grande events out to 500 m can be seen 
[3]. 

The energies of cosmic rays seen in the radio at 
KASCADE is typically a few times 10^'^ eV. In- 
clusion of KASCADE Grande provides informa- 
tion up to 10^^ eV. Below 10^^ eV the radio signal 
vanishes in the noise. 

Clearly, the Forschungszentrum Karlsruhe, where 
LOPES is located, is not an ideal location of radio 
observations due to an enormous man-made radio 
background noise (Radio Frequency Interference, 
RFI). This can be overcome somewhat in the post- 
processing through digital filtering methods. 

Also, in addition to the LOPES antennas a few ad- 
ditional log-periodic antennas ("Christmas trees") 
with new electronics have been added in order to 
develop a self- triggering algorithm [5]. Here the 
background noise is an even more severe problem. 
On the other hand, studying self-triggering under 
these conditions, will allow one to self-trigger al- 
most everywhere else in the world as well. 

Codalema 

Parallel to LOPES the CODALEMA experiment 
was set-up[4], which initially had a complemen- 
tary approach: set up particle detectors near 
an existing radio astronomy telescope at Nancay 
(France) and try self-triggering. This had the ad- 
vantage of a radio quiet site — much better than 
LOPES — and well-calibrated antennas, but the 
disadvantage of having to calibrate a new particle- 
detector array and to trigger on a yet not under- 
stood radio signal. In the latest version the CO- 
DALEMA array is now completely independent 
and employs a set of 16 wide-band active dipoles 
aligned on two 600 meter long baselines in the 
North-South and East- West directions. The radio 
array is triggered by a ground detector array of 240 
meters square containing 13 plastic scintillator sta- 
tions. The recording bandwidth is 1-200 MHz, but 
signal detection is mainly done around 50 MHz, 



where also LOPES operates. Hence, the two ex- 
periments are now quite compatible. 

LOFAR 

A next big step in radio detection of air showers 
will be the LOFAR array which was planned as 
a large radio astronomy experiment [22, 58]. In 
the summer 2007 the project had to be downsized 
due to financial shortfalls, however, it will still be 
a major step forward. According to the current 
plans, LOFAR will consist of 36 antenna fields 
("stations") in the Netherlands plus a number of 
stations across Europe (E-LOFAR: Germany, UK, 
France, Sweden, Italy, Poland, Ukraine). The first 
20 stations are expected to operate early 2009. 

Most of the antennas will be in a central concentra- 
tion ("core") of 2 km diameter, where 18 stations 
are foreseen. Each station has two sets of receiver 
systems operating from 10-90 MHz (low-band an- 
tennas, LBA, Fig. 4) and 1 10-240 MHz (high-band 
antenna tiles, HBA). The LBA field consist of 48 
dual-polarization inverted-V antennas, while the 
HBA fields consist of two sub-fields of 24 dual- 
polarization tiles. Each tile consists again of 16 
bowtie- shaped fat dipoles (Fig. 5). 

This means that in the inner 2 km there will be 
more than 800 dual-polarization low-frequency an- 
tennas and about ~ 14, 000 high-frequency anten- 
nas (> 800 tiles) that will be able to observe bright 
radio events and deliver unprecedented detailed in- 
formation about radio shower properties. LB As 
and HBAs share one receiver, so each LBA/HBA 
pair cannot observe at the same time, however, it is 
well-possible to have one half of the receivers ob- 
serve at the low-frequencies, while the other half 
observes at higher frequencies. 

For normal radio astronomical observations the ra- 
dio data from one station is combined into one 
data-stream ("digital beam-forming") to look in a 
predetermined direction. However, every antenna 
is also connected to a one second ring-buffer with 
FPGA-based processing and triggering capability. 
This allows one to trigger on the raw radio data 
stream in an intelligent way. Since 8 antennas 
share one memory board (with 4 FPGAs), there is 
the possibility to trade the number of antennas for 
triggering power or buffer length. 
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Figure 3: Left: LOPES and KASCADE Grande Layout. Right: Schematic of LOPES Electronics 





Figure 4: LOFAR Low-Band Antennas 

SKA 

As a next step in radio astronomy at the 
low frequencies the Square Kilometre Ar- 
ray (SKA) project is planned for > 2015 
(www.skatelescope.org), which will provide even 
more opportunities for radio detection of cosmic 
rays and neutrinos [21]. The SKA will employ a 
mix of receptor technologies (dipoles, tiles, small 
dishes) depending on the frequency range (70 
MHz- 10 GHz) and have a phased roll-out. The 
first phases will concentrate below 1 GHz and will 
be of high interest for radio particle detection as 
discussed here. 



Figure 5: Inside of a LOFAR High-Band Antenna 
(HBA) tile, showing one crossed bow-tie antenna. 

Auger Radio 

In the spirit of the LOPES experiment, putting ra- 
dio antennas next to existing cosmic ray experi- 
ments, it makes sense to also place radio anten- 
nas at today's largest cosmic ray array, the Pierre 
Auger observatory [66, 8]. This has been at- 
tempted recently with a few prototype radio an- 
tennas (van den Berg et al. [65]), including some 
LOPES antennas. The medium-term goal is to 
cover a 20 km^ region with self-triggering radio 
antennas. In the same region an infill aiTay and 
an upgrade to a fluorescence telescope will lower 
the energy threshold of Auger. This will nicely 
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connect to the LOPES and CODALEMA measure- 
ments in energy and allow for the first time triple 
coincidences between radio, particle detectors, and 
fluorescence and hopefully further dramatically in- 
crease the quality of the data. 

Radio in Ice and from the Moon 

Apart from the radio air shower experiments a high 
level of attention has also been attracted by the pos- 
sibility to detect showers generated in solid me- 
dia, such as ice, salt or the lunar regolith. This 
goes back to the original suggestion by Askaryan 
[6, 7]. Further theoretical progress in the 1990's 
by Alvarez-Muniz and Zas et al. [2, 69, 12] and 
successful accelerator experiments, validating the 
theory, breathed new life into this field [60, 24] . 
One idea is to use the huge detector volume of the 
moon and observe it with sensitive ground-based 
radio telescopes in search for nanosecond pulses 
which are dispersed by the Earth ionosphere. First 
such experiments were made with the Parkes ra- 
dio telescope [29] and the 64 m Kalyazin radio 
telescope [13]. An experiment (GLUE) using the 
NASA deep space network antenna at Goldstone 
received wide attention and produced interesting 
limits on the ultra-high-energy neutrinos flux [25]. 
Currently extensive experiments, LUNASKA and 
NuMoon, are progressing at the Australia Tele- 
scope Compact Array (ATCA, James et al. [38]) 
and the Westerbork Synthesis Radio Telescope 
(WSRT, Scholten et al. [62]). Later LOFAR and 
the SKA could be used to further improve the cur- 
rent upper limits to interesting levels or even detec- 
tions [61,21]. Even a detection of neutrinos using 
radio on the moon has been considered [40]. 

Instead of looking up, one can also look down on 
the Earth and use for example the large ice sheet 
of Antarctica or salt domes as detector targets. 
While salt domes [51] are currently not pushed 
very strongly due to high drilling costs, experi- 
ments involving Antarctica are flourishing. 

For quite some time already the RICE experiment 
[46] has radio antennas in the ice near the IceCube 
array [42] and a major extension of IceCube with 
radio antennas is actively discussed [42, 43]. Al- 
ternatively radio antennas have been tested at the 



Ross ice shelve in Antarctica (ARIANNA) which 
is logistically more conveniently located [9, 11]. 

An alternative approach to embedding a large num- 
ber of radio antennas in the detector volume is 
to just fly over it and use the long range capa- 
bility of radio detection. Gorham et al.[48] have 
used a military satellite (FORTE) to search for neu- 
trino induced radio pulses from the ice. Recently 
the dedicated balloon experiment ANITA [10] was 
launched for the first time to circle Antarctica and 
to detect there distant radio pulses from up-going 
neutrinos. Unfortunately the flight was cut short by 
unfortunate wind conditions, but nonetheless the 
data analysis is proceeding. 

This brief summary already shows that the number 
and breadth of radio experiments for cosmic ray 
and neutrino detection is rather large already. 

Experimental Results & Calibration 

In the following we will summarize some of 
the important experimental conclusions concern- 
ing the air shower radio properties that have been 
found recently, here mainly focused on air showers 
and based on the LOPES results. 

First of all one has to realize that the simplicity of 
the antenna comes at the cost of more complicated 
calibration. The sensitivity of a dipole depends on 
frequency, direction, and polarization. For LOPES 
the absolute calibration has been performed as part 
of the PhD thesis of S. Nehls [54] using an ele- 
vated calibrated reference antenna. On the other 
hand, given the simple structure of the antenna, it 
is also possible to calculate the expected beam pat- 
tern on the sky using standard antenna simulation 
packages. 

An example is shown on Fig. 6, which shows that 
the beam shape is elongated and even not peak- 
ing towards the zenith for frequencies above the 
resonance frequency of the LOPES dipoles (^60 
MHz). The elongated structure is related to the ori- 
entation of the dipoles (EW) and would be rotated 
by 90° for the other (NS) polarization. In turn this 
also means that the crossed-dipole, if uncalibrated, 
will always produce highly polarized signals. 
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Figure 6: Gain pattern of a LOPES antenna at 50 
MHz (top) and 70 MHz (bottom) showing the en- 
tire sky. The color scale goes from a gain of 
(black)to 5 (white). 

Energy Calibration 

In addition to the antenna dependencies, the radio 
emission will also depend on the shower geome- 
try. The main factors that have been identified are 
the particle energy, Ep, the angle between shower 
axis and geomagnetic field ("geomagnetic angle"), 
a, the zenith angle, 9', and the distance from the 
shower core, r. For an east-west polarized an- 
tenna one finds that the radio emission for showers 



at the same energy and distance is proportional to 
1 — cos a (Fig. 7a), the signal drops exponentially 
with radius (Fig. 7b), and increases linearly with 
primary particle energy (Fig. 7c). 

Altogether this has been nicely parametrized in 
Horneffer's formula [3 1 ] (at the moment valid only 
for the EW polarization): 

Cost = (11 ±1.) ((1.16 ±0.025) -COS a) COS 61 
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We note that the exponential decay of the radio sig- 
nal is also seen by the CODALEMA experiment 
[4]. 

This prescription can now be inverted to predict 
the energy of the incoming particle. Comparison 
between the energy predicted from radio with the 
energy estimated from KASCADE Grande, shows 
a scatter of 27% between the two methods for 
Ep > 10^^ eV. This is very encouraging, given that 
shower-to-shower fluctuations in the KASCADE 
Grande estimate alone should produce a 25% scat- 
ter. Hence, the scatter in the radio measurements 
should be much less. 

This would support the claims from Monte Carlo 
simulations [35] that radio is a good tracer of the 
energy. The main reason why one suspects lower 
scatter in the radio measurements with respect to 
particle detection on the ground is the fact that ra- 
dio emission in not absorbed in the atmosphere. 
Hence, radiation from every particle is visible on 
the ground — it is in that sense a bolometric mea- 
surement. Variations in the location of the shower 
maximum will be less dramatic compared to mea- 
surements of particles on the ground which are just 
a fractional tail of a quickly declining function, 
whose values are quite sensitive to X^ax- 



Spectrum 

One topic that had been difficult to tackle in the 
past has been the spectral shape of the radio sig- 
nal, i.e. how much power is emitted at which fre- 
quency? Historic experiments were relatively nar- 
row band and non-simultaneous data had to be 
combined. Modern broad-band receivers allow one 
to study the instantaneous spectral index, but re- 
quire careful bandpass calibration. 
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Figure 7: Calibration results of LOPES showing 
the normalized radio voltage vs. air shower param- 
eters where the other parameters have been divided 
out. Each point represents the average and spread 
of all events in that bin. Panels a-c are from top to 
bottom. 

Nigl et al. [56] employed two methods to get the 
spectral shape (Fig. 8): Fourier transform of the 
(not squared) electric field around the pulse posi- 
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Figure 8: Example spectrum of a LOPES event for 
two different methods of beam-forming. The mea- 
sured spectral slope is e oc v^^. The lines rep- 
resent two spectra from Monte Carlo simulations. 
Also shown is the noise spectrum that has been cor- 
rected. 



tion and measurements of the pulse heights after 
applying narrow-band digital filters to the signal. 
Both methods give consistent results. 

The spectra can be represented by a power-law 
function or an exponential decay. For the narrow 
frequency range of LOPES, extending only a factor 
of two, we cannot distinguish between the two pre- 
scriptions. For a power-law function the average 
spectral index is i^^^^" •^. This would mean that 
the power of the signal falls of with i'^^. This is 
consistent with the simple expectations of coherent 
geosynchrotron [20] and only slightly steeper than 
the Monte Carlo simulations suggest [34]. Also, 
the Codalema experiment finds powerlaw spectra 
with spectral indices in the range -1.5 to 0. LOPES 
sees spectral indices in the range -1.5 to -0.4, so 
there is some agreement, but more detailed inves- 
tigations have to be performed in the future. 

Direction & Imaging 

The next question then is, how well can we localize 
the radio emission? This has become of particular 
importance given the finding of anisotropies and 
correlations between cosmic ray arrival directions 
and nearby extragalactic objects [8]. 
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Figure 9: Image of a bright radio flash at the time 
of maximum. The dark dots mark the fitted LOPES 
position and the direction found by KASCADE. 
The size of the ellipse is determined by the ex- 
pected image resolution. 

Using radio astronomical imaging techniques, we 
can actually image the radio flash from the air 
shower. For LOPES L. Bahren has developed a 
special tool ("skymapper") which can actually do 
this on the tens of nanosecond (i.e., the sampling 
rate) level (see Fig. 9) and in three dimensions 
(Fig. 10). 

This issue of the positional accuracy has then 
been further investigated by Nigl et al. [55] us- 
ing LOPES data. Conventional interferometry is 
very sensitive to positional changes. For point 
sources one expects an angular error Aamm — 
±igj|p--^ in the azimuthal direction, where SNR 
is the signal-to-noise ratio, D the separation of the 
antennas, and A the observing wavelength. For 
high SNR images the positional accuracy for point 
sources is always better than the image resolution. 
Hence, for an SNR of 10 for an antenna separa- 
tion of 100 m and observing frequencies around 60 
MHz, as in LOPES, one expects an error of only 
0.15°, while the point-spread function of the inter- 
ferometer has a much larger width of about 2-3° in 
azimuth. 



Figure 10: Three dimensional view of the radio 
emission from a cosmic ray air shower detected 
with LOPES. 3D sidelobes have been "cleaned" by 
only displaying high-brightness regions. The peak 
at 3 km height appears to be real, but the extended 
structure along the axis is strongly affected by re- 
maining sidelobes. 



Comparisons of the shower direction (assuming 
a fixed shower core) between KASCADE and 
LOPES, actually shows an average offset of 1.3° 
(Fig. 11), which does decrease with increasing sig- 
nal level (and increasing SNR). This is not bad 
compared to the imaging resolution of LOPES, but 
we would have expected better results still. So, 
what is the dominating source of error? 

Some insight can be obtained from Fig. 12 which 
shows that the location of the radio centroid on the 
sky is also a function of the distance or radius of 
curvature. One has to remember that the shower 
maximum of the showers that LOPES sees is just a 
few km high. This is still in the near field of the in- 
terferometer. This means that radio waves emitted 
in the shower maximum will not appear as a plane 
wave, but will show a curvature with a radius cor- 
responding to the distance from the observer In the 
real world, the wave front will be even more com- 
plicated, since the emission is not constrained to 
a small region but extends along the shower axis. 
In principle one has then a wavefront that is the 
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Figure 11: Offsets on the sky between the di- 
rection of the air shower axis as measured with 
KASCADE and LOPES vs the radio pulse ampli- 
tude. Events with circles were taken during thun- 
derstorm conditions. 
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Figure 12; Cross section along the elevation axis in 
an image of a cosmic ray radio flash for different 
radii of curvature. The maximum field strength is 
found around a radius of 3km. 



superposition of many spherical waves emitted at 
different positions. 

In a proper 3D imaging process one would try to 
deconvolve the data and put together a 3D image 
cube. We are not yet able to do this. So, all we can 
do at present is to try different radii of curvature 
and search for the maximum in the emission. This 
is what Fig. 12 shows: different cross sections of a 
radio image focused at different distances in steps 
of 250m. The maximum is found at a distance of 
about 3 km. This is 30 times farther than the typi- 
cal baselines on the ground and only possible due 
to the good SNR. What is clear from the figure is 
that not only the emission level changes but also 
the position of the maximum. The problem is that 
for small radii of curvature a small change in ra- 
dius is similar to an inclination of a plane wave. 
Inclining the (virtual) receptor plane implies a po- 
sitional shift on the sky. As seen in this example, 
the shift can be up to 3°. Hence, an error in the ra- 
dius of curvature determination or - perhaps more 
important - a non-spherical wavefront will propa- 
gate into an positional error! 

This requires that radio shower parameters are re- 
ally derived from a 4 (or 6) dimensional data cube, 
consisting of time, 3 spatial coordinates, and po- 
tentially 2 shower core location parameters [3]. 



So, any experiment will improve its spatial accu- 
racy not only with greater baselines but also with 
increasing the number of antennas. In addition we 
need to understand the exact geometry of the radio 
shower front. Here, further simulations and the de- 
tailed observations with the many antennas of LO- 
FAR should clarify that issue. This may make fur- 
ther dramatic improvements in the astrometry of 
radio air showers possible. 

Electric Fields and Lightning 

One other important factor that has been looked at 
with some worries, is the influence of the atmo- 
spheric electric field on the radio emission. While 
the Earth magnetic field is very stable, the electric 
field can change significantly from 1-10 V/cm dur- 
ing fair weather, to 100 V/cm in heavy rain clouds 
(Nimbostratus), and up to 1000 V/cm in severe 
thunderstorms. If the radio emission is affected by 
the electric field one would not be able to inter- 
pret the radio signal quantitatively, since measur- 
ing precisely the instantaneous electric field struc- 
ture is almost impossible. 

To get a first idea of the importance of the electric 
field, we could simply look at the Lorentz force, 
F — e{E + v/c X B) in cgs, which is driving the 
geosynchrotron emission. For a relativistic parti- 
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Figure 13: Radio pulse height (normalized by en- 
ergy) versus the geomagnetic angle. Blue dots are 
taken during clear weather, purple during nimbo- 
stratus, and red points during thunderstorm condi- 
tions. 

cle (v/c ~ 1) the electric field will then dominate 
if £; > B ~ 150Vcm-i(S/0.5G). Hence, from 
this simple approximation one would expect only 
for severe weather a modification of the radio sig- 
nal. 

Buitink et al. [16] investigated radio pulse heights 
with LOPES under different weather conditions. 
They selected time slots where the local weather 
station recorded clear weather, heavy rain, or 
thunderstorms and compared the three (Fig. 13). 
Within the errors, the radio pulses followed the 
previously found correlation with the geomagnetic 
angle (if normalized to the same energy and ra- 
dius). However, significant outliers are found in 
the thunderstorm data set - and only there. The 
amplification of the radio signal can be a factor ten. 

This seems to confirm the simple estimate we made 
above and the fact that the Lorentz force is dom- 
inating. Further verification comes from Monte 
Carlo simulations. Buitink et al. [17] have now in- 
cluded electric fields in the CORSICA and REAS2 
codes, allowing one to model the E-field influ- 
ence in detail. The simulations have calculated ra- 
dio pulses for different values of the electric field. 
Again a significant amplification is seen as soon as 
flie E-field reaches 1000 V/cm, while at 100 V/cm 
the radio emission shows only little differences. 



The CORSIKA simulations also show a modifi- 
cation of the electron/positron energy distribution 
and the shower structure. Many pairs are deflected 
significantly from the shower axis, which might be 
detectable with particle detectors. 

Finally, we note that besides the radio pulse height 
also other parameters are impacted by thunder- 
storm electric fields. For example, in the investiga- 
tion of the emitted radio spectrum one bright radio 
event stood out with a much steeper radio spectral 
index. It was found to be a thunderstorm event. 

Moreover, also in the positional offsets between 
LOPES and KASCADE thunderstorm events stand 
out. In Fig. 1 1 all outliers are events measured dur- 
ing thunderstorms. Whether this is due to an actual 
deflection or an asymmetry in the radio emission is 
not yet clear 

In summary, we can state that air showers passing 
through the strong fields of thunderstorm clouds 
are brighter, further offset from the shower core 
measured on the ground, and have a steeper radio 
spectrum. For the measurement of cosmic rays this 
means that radio — at least at frequencies above 
40 MHz — remains a rehable technique as along 
as data taken during thunderstorm is discarded. 

On the other hand, the current results strongly sug- 
gest that in high E-fields not just the radio emis- 
sion is altered, but the entire shower (at least the 
electronic part). This point may warrant further 
investigation for its own sake. Moreover, it has 
been speculated that cosmic ray air showers could 
play a role in initiating lightning through a run- 
away breakdown effect [28]. 

Radio methods could help to investigate this con- 
nection experimentally, since both - the lightning 
strike and the air shower - would be detectable 
by the same instrument. Also, further Monte 
Carlo simulations will investigate whether there 
is enough energy gain through the electric field 
or ionization through the air shower to actually 
start the runaway breakdown process. The LOFAR 
project will try to address some of these issues. 

Self-triggering 

Overall, the current results have provided a very 
comprehensive picture of the radio properties. 
However, whether the radio detection will mature 
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into a standard technique depends on whether it 
is possible to actually trigger on the radio sig- 
nal. This has been attempted but the final break- 
through still stands out. Recent attempts were 
made for example by the CODALEMA experi- 
ment [57]. Within LOPES, LOPESS"^^^ [5] has 
been designed specifically to investigate this issue 
in more detail. 

We do now understand where the challenges lie. 
In many cases, disturbing radio interference has 
actually been generated by the devices and elec- 
tronics of the cosmic ray experiments themselves. 
So, designing radio quiet electronics, power sup- 
plies and data communication is an important first 
step. Also, radio contains an enormous amount of 
information (as witnessed by every FM receiver) 
and many processes of the modern world produce 
wanted and unwanted radio signals. So, filtering 
out the correct information is more complicated 
then just looking for a peak in the electric field. 

In a broad-band receiver the biggest contribution of 
the basic noise level typically comes from narrow- 
band RFI transmitters. Hence, a digital filter to 
cut these signals out - which is currently only em- 
ployed in the post-processing, is a crucial step in 
the triggering electronics. Moreover, the charac- 
teristics of the pulse itself need to be considered 
as well. For a human eye it is quite simple to dis- 
tinguish a cosmic ray pulse from those generated 
by a passing 1970 Chevrolet. Hence pulse shape 
parameters need to be used in the triggering. This 
will in any case require a bit more intelligence on 
the trigger board, than with conventional experi- 
ments. Such a pulse-shape parameter search to im- 
plement self-triggering in hardware is currently un- 
der way at LOPES ^^^'^ showing some interesting 
progress recently. These techniques will eventu- 
ally be tested at Auger and one can be hopeful that 
this last and crucial step will be achieved in the not 
too distant future. 

Radio Pulses from the Moon 

The main focus of this article was on the detec- 
tion of radio emission of air showers. However, we 
want to end with a few comments on the prospects 
for radio emission from ultra-high energy cosmic 
rays with upcoming radio telescopes, in particular 



with LOFAR. Here the atmospheric detection will 
be implemented in the project — with the "tran- 
sient buffer board" being at the heart of this new 
technique. It turns out that this buffer-board may 
also improve the detection of cosmic rays hitting 
the moon. 

Scholten et al. [61] have shown that the 100-200 
MHz range is ideal for detecting cosmic rays above 
lO^o eV hitting the lunar surface. In LOFAR any 
such event could be detected in a beam formed 
towards the moon. This detection could be used 
to trigger the buffer boards and download the raw 
data from all LOFAR antennas. With the raw data 
of the individual antennas the exact nature and ori- 
gin of the pulse could be determined much more 
precisely, if one uses some of the offline process- 
ing steps known from the air shower detection. 

The question is how sensitive is this technique? 
The originally planned LOFAR would have been 
very favorable for this [61], however, due to the 
downsize of LOFAR the sensitivity has decreased. 
K. Singh has now recalculated the expected sensi- 
tivity with the latest available station layout. The 
result is shown in Fig. 14. With the lower sensitiv- 
ity, the minimum energy that can be detected has 
moved up above 10^^ eV, which will decrease the 
expected count rate. 

For comparison we also show the latest Auger 
spectrum with an extrapolation of the power law 
with three spectral indices through the GZK cut- 
off. For a powerlaw of E^^ LOFAR could in prin- 
ciple reach this extrapolation with a total observ- 
ing time of 90 days, if a large number of tied-array 
beams can be formed - a mode that still needs to 
be tested. 

While strong evidence for a GZK cut-off above 
j^q19.6 gY jj^g been seen by Auger, it cannot be 
excluded that there will be some recovery of the 
spectrum from local sources or neutrinos. Hence, 
it is still worth looking in this regime. Radio will 
probably be the only technique that can deliver at 
least very meaningful upper limits in this energy 
range. 

Future expansions of LOFAR, such as the SKA, 
can improve these limits even more, actually reach- 
ing down to the GZK cut-off with the clear expec- 
tation of actual detections. This will then proba- 
bly provide the ultimate detection experiment for 
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Figure 14: Expected LOFAR sensitivities for detecting ultra-high energy cosmic rays hitting the moon. "Full 
stations" here relate to the originally planned 96 tiles. 36 quarter stations are the currently planned layout 
for the central LOFAR core. "Coherent" and "incoherent" relate to the two methods of combining the data 
from different stations. 



cosmic ray particles at the highest energies we can 
ever measure. 

Conclusions 

The radio detection technique for cosmic particles 
has seen quite some ups and downs in the last 
decades. Pronounced dead in the 1970's it has now 
risen from the ashes. But is it here to stay? The 
chances are good at least. 

First of all, we have made major progress in un- 
derstanding the emission mechanism. For radio 
in solid media (Cherenkov emission) codes have 
been developed and accelerator experiments have 
been performed, giving trust in the reality of the ef- 
fect. For radio emission from air showers (geosyn- 
chrotron) good Monte Carlo codes and solid ex- 
perimental verification are now available and more 
and more details are being worked in. 

The LOPES experiment, and in some areas also 
CODALEMA, has given us already detailed infor- 
mation about the radio air shower properties and 
performed a very useful cross-calibration between 
particle and radio detectors. Moreover, major ex- 



periments have embraced the technique. The large 
LOFAR radio telescope has the cosmic ray detec- 
tion built in, the Auger collaboration is testing it, 
ANITA has flown, and also the IceCube collabora- 
tion is seriously preparing radio experiments. 

A few issues still need to be solved: how does an 
optimal trigger system look like for radio anten- 
nas? What is the optimal layout for a large radio 
array? Nonetheless, there is a good chance that ra- 
dio detection will become common place over the 
next few years. 

It will be interesting to see in which direction this 
will develop. The simulations for air showers indi- 
cate that the addition of radio may increase energy 
resolution and directional accuracy. Also compo- 
sition information (through Xmax) seems to be en- 
coded in the radio signal. After the breakthrough 
of the hybrid technique with Auger, maybe we will 
see "tri-brid" detectors in the future. More and 
better data is almost always better in physics. If 
cosmic ray air shower arrays are to continue in the 
next decades, they will likely include radio anten- 
nas. For the highest energy events, the new gen- 
eration of low-frequency radio telescopes provides 
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hope that even particles above lO^'^ eV could in 
principle be detected in the future — if they exist. 
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